Intrinsic Josephson effect and nonequilibrium soliton structures in two-gap 

superconductors 
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We predict a new dynamic state in current-carrying superconductors with multicomponent order 
parameter. If the current density J exceeds a critical value Jt, an interband breakdown caused by 
charge imbalance of nonequiUbrium quasiparticles occurs. For J > Jt, the electric field penetrating 
from current leads gives rise to various static and dynamic soliton phase textures, and voltage 
oscillations similar to the nonstationary Josephson effect. We propose experiments to observe these 
effects which would probe the multicomponent nature of the superconducting order parameter. 

PACS numbers: PACS numbers: 74.20. De, 74. 20. Hi, 74.60.-w 



There are experimental evidences that several super- 
conductors, including MgB2 0, ^, NbSe2 @, the heavy- 
fermion UPt^ 0, the organic {TMTSF)2X g and 
K — BEDT compounds may have a multicomponent 
order parameter with internal degrees of freedom ||^, ^ . 
For instance, for two weakly coupled s-wave order param- 
eters -01 = Aie**i and ■02 = A2e'^2 corresponding to two 
different parts of the Fermi surface (like in MgB2), the 
internal degree of freedom is the interband phase differ- 
ence 0{r,t) = 6i — 02- In this case, besides the usual 
phase-locked states {9 ~ or tt), intrinsic phase tex- 
tures 0{r,t) can occur. The soft modes associated with 
the fluctuations of 0{r,t) are nearly decoupled from the 
gap fluctuations and behave like the Anderson plasmons 
in Josephson junctions Q. These modes may manifest 
themselves as resonances in the ac Josephson effect 
or static 27r phase kinks in 0{x) [ p^ . 

In this Letter we show that, although static metastable 
textures 9{r) do not manifest themselves in the equilib- 
rium magnetic response, the situation radically changes 
for nonequilibrium current states where the charge im- 
balance near normal leads gives rise to dynamic phase 
slip structures 0(y,t) propagating into a superconductor 
(Fig. 1). We predict a new dynamic state above the 
critical current density J > Jt which marks the onset 
of the current-induced breakdown of the superconduct- 
ing state. For the multicomponent 0, such breakdown 
is a two-stage process. First, at J — Jt, an interband 
phase breakdown occurs, resulting in spontaneous phase 
solitons in d{x, t) and ac voltage oscillations. At higher 
currents, J = Jd, both gaps Ai^2 are suppressed by the 
pairbreaking effects. For weak interband coupling, Jt is 
much smaller than the depairing current density Jd- 

We obtain the equations of motion for 9 and the elec- 
tric field E near the critical temperature T k Tc, using 
the time-dependent Ginzburg-Landau (TDGL) equations 
Ijnj generalized to a two-gap superconductor: 



where /i runs from 1 to 2, is the electric potential, 0o is 
the flux quantum, c is the speed of light, F^ are damping 
constants, and the free energy F = J (Pr{fi + f2 + fm + 
fint) contains the magnetic part /,„ = |V x Ap/87r, the 
GL intraband part and the interband interaction fint 
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A 0. 



fint = 7(V'102 + 01 ^^2) = 27A1A2 cos 6*. 



, (2) 
(3) 



Here A is the vector potential, 7 is a coupling constant, 
the asterisk means complex conjugation. The qualita- 
tive results of this work remain valid for any other pe- 
riodic dependencies of fmtiP) on 6, for example, fint — 
2A^A2(7i cos 26* + 72), suggested for heavy-fermion su- 
perconductors Q. We focus here on two-gap supercon- 
ductors for which two cases are possible: (i) The Cooper 
instability occurs at the same temperature in both bands, 
thus a^{T) oc T — Tc, 7 oc T — Tc, while g^, and (3^ 
are constants, (ii) The critical temperature Td in the 
band 1 is higher than Tc2, in the band 2, so 0i appears 
spontaneously below Tc, while -02 at T^i < T < Tc is in- 
duced by the interband interaction, because ai (x T—Td, 
a2 oc T — Tc2, and 7, and are constants In 
this paper we consider the first scenario, assuming that 
7 ^ 0:1,2, which is hkely to occur in MgB2 ||l|, §]. 
Eqs. (|l|) yield the following gap equations 

F^A^ = .g^V^A^ - d^A^ - /3^A3 - 7A^cos6l, (4) 

where /I = 2 if /i = 1, and JI = 1 ii /i = 2, 9 = 61 — 62, 
5^ = + (27r(3p/0o)^5p: and = A - 0oV6'^/27r. 
The imaginary part of Eq. (|^) gives 

- " ^"%)A^=-^V(A^Q,)±7AiA2sin0(5) 



Po ' <Po 

V X V X A = (47r/c)(crE + J^) 



(6) 



r^(9t'0^ - 27rd(p/0o) = -SF/Sip^ 



(1) 



with the plus sign corresponding to ^ — 1. The current 
density J = crE -I- has the ohmic component J„ = crE 
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proportional to the electric field E = — V(/5 — A/c, and 
the normal conductivity a. The supercurrent 



-8^2^(51 A2Qi+.g2A2Q2)/02 



(7) 



is a sum of independent intraband contributions. The 
boundary conditions, + 27rA„/(/)o)'0/x = iV'/x/^A" 
tween a superconductor and a normal metal imply zero 
perpendicular supercurrents for both -01 and ip2 ■ Exclud- 
ing Qi,2 = -((/'oJ. ± 47rc524A|iV6l)<?!)o/[87r2c((?iA2 + 
52 from Eqs. (||) and (|^), we obtain the energy den- 
sity f = fe + fe + E^KK + ffM(VA^)2 + /3^a4/2], 
where the magnetic and phase dependent parts, fm and 
fe turn out to be decoupled from each other: 



jj2 



fe 



27A1A2 cosf 



(8) 
(9) 



For constant gaps A^ , variation of Eq. (||) with respect 
to H yields the London equation A^V^H — H with the 
magnetic penetration depth — 0Q/327r'^(5iA^ -I-52A2). 
If 7 ^ ai^2, the gaps A^ are decoupled from 9, so varia- 
tion of Eq. (^) with respect to 9 yields the sine-Gordon 
equation LgW^O = ± s'mO which has a single-soliton solu- 
tion ||l^ , or periodic solutions analogous to the Josephson 
vortices in a magnetic field [T^ . However, the interband 
phase difference 9 is decoupled from static magnetic fields, 
while the Gibbs free energy of a Josephson contact con- 
tains the interaction term oc (HV^j). Thus, equilibrium 
nonuniform solutions 9{x) are energetically unfavorable 
as compared to the uniform phase-locked state 9 — 
for 7 < or 6* = TT for 7 > 0. By contrast, various 
dynamic or quenched phase textures can be generated 
during current-induced interband breakdown. 

We consider the weak coupling limit 7^0;^, for which 
the uniform gaps Ai^2 are unaffected by 9{r,t), so Eqs. 
®"® S^^^ ^ time-dependent London equation with the 
account of dissipation and charge imbalance effects . 
Expressing Q1.2 in Eqs. in terms of J and V0, and 
then subtracting the equations for 9i and 92 from each 
other, we obtain the following equation for 9 



Tg9 = igV^6l ± sin6' -f a^divJ^ 



(10) 



where the relaxation time rg, the decay length Lg, and 
the charge coupling parameter ag are given by 



Tg = rir2AiA2/|7l(riA2 + r2A^ 

:5ig2AiA2/|7|(5iA?+g2A^ 



Ll 



ag 



/)oAiA2(rig2 - r2.gi) 



4^c|7|(giA2+52A2)(riAf +r2A2) 



(11) 
(12) 

(13) 



The signs in Eq. (^^ correspond to the sign of 7. As seen 
from Eq. (p^, the phase mode is unaffected by any dis- 
tribution of bulk supercurrents, unless there is a nonequi- 
librium charge imbalance divJs = — crdivE caused by the 



electric field penetrating from normal current leads. The 
excess charge density provides a driving term in the sine- 
Gordon equation for 9, where the coupling constant ag is 
nonzero if giT2 ^ 52rij which implies different electron 
diffusivities in the bands 1 and 2. 

To obtain the equation for E, we add Eqs. (||) for 9\ 
and 02) then take the gradient of the sum and express 
Qi,2 in terms of J and V6'. This yields 



TeE + E - LggraddivE + a^\I9 = j/cr. 



(14) 



where J(i) is the driving current density, is the electric 
field penetration depth, Tg is the charging time constant, 
and the coupling term a^9 describes an electric field 
caused by moving phase textures: 
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f7</.;^/87r^c^(riA2 4 
Ue = 2|7|AiA2Q!e 



92AI), 

^2Al), 



(15) 
(16) 
(17) 



We first use Eqs. ( |10| ) and (^_4|) to calculate the phase 
textures in a microbridge of length 2a (Fig. la) for which 
Eq. ( p^ ) at 7 < takes the form 



Ll9" 



sine + agj' = 



(18) 



Here 9'{x) and Js{x) vanishes at the bridge edges, x = 
±a, where E{x,t) — i?o cosh(a;/Le)/ cosh(a/Le) equals 
the value = J j a provided by an external power 
source, 9' (0) — E'{0) = 0, and the prime denotes dif- 
ferentiation over X. Now we obtain the condition un- 
der which Eq. ( |l8| ) has a stable solution 9{x,9a) lo- 
calized at the edge of a long (a ^ Lg) bridge, where 
9a — 9{a). We first consider the limit Lg ^ L^, for which 
Jg{x) = —aE'{x) is essential only in a narrow region 
a — Le < X < a. Multiplying Eq. ( |l8| ) by 9' and integrat- 
ing from to a — Lg, we obtain 9'{a — L^) — — 2sin6'a/2. 
Here 9'{a—Le) — —agJ/Lg is obtained by integrating Eq. 
( p8| ) over a — L^ < x < a, where 9{x — Le) ~ 9a, and sin9 
can be neglected. Excluding 9' {a — Lg), we arrive at the 
equation agJ — 2Le| sin6'a/2|, which has solutions for 9a 
only below the critical current density Jt = 2Lg /ag : 



_ STTciT^Al + r2Ai) V|7|gig2(giAj + ff2A^) ^^^^ 
(f>o^/AiA2\giT2 - 52ri| 

In the opposite limit Lg <C Le, the Laplacian 
in Eq. ( |lo| ) can be neglected. Then the steady- 
state solutions sm9{x) = —aaeE'{x), and E{x) = 
i?o cosh(a;/Le)/ cosh(a/Le), exist only if aaeE'{a) < 1, 
or J < Jt = Lg/ofe tanh(a/Le), where 



Jt = 



[a(riA2+r2A2)]l/2(g^A2+52A2) 



(20) 



V2AiA2|.gir2 - .g2ri| tanh(a/Le) 

If 7(T) and ^^(T) linearly pass through zero at Tc, Eqs 
M) and M) yield Jt cx (T^ - T)3/2. 
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For J > Jt, the charge-induced interband breakdown 
gives rise to a striking dynamic phase texture in which 
sohtons periodically appear near the current leads and 
then propagate to the bulk. The actual soliton dynamics 
depends on the particular boundary conditions, so we 
performed numerical simulations of Eqs. (^-(^ for 
the different geometries shown in Fig. 1 . We consider the 
limit aeag <C min(reLg, TgLg) of weak coupling between 
9 and E, for which the driving term divE in Eq. ( p^ ) 
is mostly determined by the static E{x) near the current 
lead, while the last term in Eq. gives only a small 
ac correction to E of the second order in . 

Fig. 2 shows the soliton formation near the current 
lead in the thin microbridge (Fig la) as J was instanta- 
neously turned on from zero to a value below Jt- Such 
current step produces a phase soliton localized near the 
edge, while the bulk of the bridge remains in the phase- 
locked state with either 6* = or tt, depending on the sign 
of 7. For J < Jt, this behavior is characteristic of all ge- 
ometries shown in Fig. 1, whereas the dynamic phase 
texture for J > Jt can be very different. 

We start with the bridge geometry (Fig. la), for which 
E(x,t) and 9(x,t) are even and odd functions of x, re- 
spectively, so the boundary conditions are: E{±a, t) — 
Eq, E'{0,t) = 0, 61(0) = 0, and B' (±a,t) = 0. The lat- 
ter condition ensures that supercurrents in both bands 
vanish at the normal electrodes, where J ~ uEq. In 
this case a soliton first appears at the bridge edge, but 
for J > Jt, it is pushed to the bulk by the strong gra- 
dient of E{x). Then the next soliton forms near the 
edge and the process repeats periodically, resulting in 
the propagation of the soliton chain into the bulk, as 
shown in Fig. 3. After the first soliton reaches the 
center of the bridge, it stops (because 9{Q) = 0), while 
new solitons keep entering the bridge from the current 
lead. During this soliton pileup, the mean slope 9'{t) 
increases, reaching a critical value 9' at which the soli- 
ton generation at the edge stops and a static phase tex- 
ture forms. For J 3> Jc, the maximum soliton density 
n = can be estimated from the static Eq. ( |l0| ) 

in which the rapidly oscillating sinB can be neglected, 
and E{x) — J smh{x/ Le)/(T smh{a/ L^)- This yields 
9'{x) = aoJLe[cosh.{a/Le) — cosh(a;/Le)]/Lg sinh(a/ie), 
giving for the long bridge (a L^): 

9', - JLe<^o(ri.g2 - T2gi)/A^cgig2{Ti^l + FaA^) (21) 

During the formation of the soliton chain, t < 
tc ~ Tga9'^/2TT, voltage oscillations are generated on 
the bridge. A similar transient behavior occurs at 
the point contact (Fig. lb), for which divE = 
21 exp{—r / Le)rQ/TTLe{rQ -\- Lf.)r, where / is the current 
through the semi-spherical contact of radius rg . In this 
case the increase of I above Ic causes propagation of con- 
centric soliton shells until the static structure forms as 
the critical phase gradient ~ 9'^ is reached. 



A very different soliton dynamics occurs in the 4- 
terminal geometry (Fig. Ic), for which the supercur- 
rents make 90° turns around the central stagnation point 
(x = 0) where V9 = 0. As a result, E{x) and 9{x,t) be- 
come odd and even functions of x, respectively, which 
radically changes the dynamics of the phase textures as 
compared to the bridge, since now the solitons do not 
stop as they reach the center, but keep moving until they 
reach the opposite current lead. This is due to the fact 
that the driving charge density divE does not change sign 
all the way along the horizontal leg of the cross in Fig. Ic. 
The corresponding total extra charge (which plays the 
same role here as the transport current in a Josephson 
junction) is exactly compensated by the opposite charge 
distributed along the vertical leg. Thus, the mean phase 
gradient 9' remains constant, while the phase slippage at 
the normal leads results in a continuous soliton shuttle 
and periodic voltage oscillations between the current lead 
and the center as shown in Fig. 4. For ^ Lg, the ac 
voltage Vui{t) = ^^^-^^Vmcoslmujt + tpm) can be esti- 
mated by integrating Eq. (Q) from to a, which yields 
TeK; + = ae[9ait) - 9(0, t)]. Here the values 9ait) ~ 
9{0,t) obey the RSJ equation TgBa + sinBa = f3, for which 
2 j: Tn=iiP - y/W^r' cos mcjt] , 
Tg''^\/l3^ - 1 pll. As a result, we 



Te9a = V^23T[i 
P = J/Jt, and u! 



obtain tan tpm = mujTg, and 



Te 



(22) 



VrniP) is maximum ai f3 — P„i, where Pm — 1-1- 
0.5(r6//mre)^/'^ for Tg ^ t^. Thus, a two-gap supercon- 
ducting bridge between normal banks exhibits the behav- 
ior of an overdamped Josephson junction. If the widths 
of the superconducting legs in Fig. Ic are different, either 
a unidirectional and mutually orthogonal soliton motion 
may occur, if J > Jt in one or both legs, respectively. 

The soft interband phase mode could also manifest it- 
self in rf absorption at frequencies below the small gap 
A2. Interaction of the 6'-mode with the rf electric field 
depends on the polarization of E: if E(t) is parallel to 
the sample surface, then divE = 0, so the phase mode 
is not excited by the rf field. If, however the rf field has 
a component perpendicular to the sample surface, then 
the phase mode contributes to the rf impedance. 

In conclusion, we predict an interband breakdown in- 
duced by nonequilibrium quasiparticles in two-gap super- 
conductors. It results in spontaneous generation of static 
and/or dynamic phase textures, causing the Josephson 
voltage oscillations. The observation of these oscillations 
would unambiguously indicate the multicomponent na- 
ture of the superconducting order parameter. 

This work was supported by the NSF MRSEC (DMR 
9214707), AFOSR MURI (F49620-01-1-0464) (AG), and 
by the US DOE Office of Science under contract No. 
W31-109-ENG-38 (AG and VV). 
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FIG. 1: Geometries in which the interband phase breakdown 
could occur. Here N labels normal electrodes, gray domains 
show phase solitons moving along thin arrows, and block ar- 
rows indicate current directions. Static phase textures form 
in microbridges (a) and point contacts (b), while in the four- 
terminal geometry (c) the solitons continuously move from 
one current lead to another. 
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